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Abstract 

We investigate the entropy production within dissipative hydrodynamics in 
the Israel-Stewart (IS) and Navier-Stokes theory (NS) for relativistic heavy 
ion physics applications. In particular we focus on the initial condition in 
a 0+lD Bjorken scenario, appropriate for the early longitudinal expansion 
stage of the collision. Going beyond the standard simplification of a massless 
ideal gas we consider a realistic equation of state consistently derived within 
a virial expansion. The EoS used is well in line with recent three-flavor 
QCD lattice data for the pressure, speed of sound, and interaction measure 
at nonzero temperature and vanishing chemical potential (/Xq = 0). The 
shear viscosity has been consistently calculated within this formalism using 
a kinetic approach in the ultra-relativistic regime with an explicit and sys- 
tematic evaluation of the transport cross section as function of temperature. 
We investigate the influence of the viscosity and the initial condition, i.e. 
formation time, initial temperature, and pressure anisotropy for the entropy 
production at RHIC at a/snn = 130 GeV. We find that the interplay between 
effects of the viscosity and of the realistic EoS can not be neglected in the 
reconstruction of the initial state from experimental data. Therefore, from 
the experimental findings it is very hard to derive unambiguous information 
about the initial conditions and/or the evolution of the system. 
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1. Introduction 

Understanding the rich phase structure of quantum chromodynamics for 
the density-temperature plane is a challenge for theoretical as well as expe- 
rimental high energy physics. Lattice Monte Carlo simulations have revealed 
several exciting results over the past decade [1, 2, 3, 4, 5, 6, 7, 8]. At high 
densities, where the lattice calculations can be performed, effective models 
of QCD are needed for a theoretical description of this rich phase structure 
[9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22]. From the experimental 
point of view, heavy-ion collisions are the tool for such investigations. Re- 
cent observations at the relativistic heavy-ion collider (RHIC) at Brookhaven 
National Laboratory (BNL) indicated that the quark-gluon plasma (QGP) 
created in ultra-relativistic Au -|- Au collisions is a strongly interacting sys- 
tem. Due to this experimental evidence the QGP can not be satisfactory 
described by the Stefan-Boltzmann (SB) limit for relativistic noninteracting 
massless particles, but a realistic equation of state (EoS) has to be used. 
Recently, we have systematically derived such an EoS within a virial ex- 
pansion [22]. The recent three- flavor lattice QCD data are described very 
well for the main thermodynamics quantities, e.g., pressure, entropy density, 
speed of sound, and interaction measure. Regarding heavy ion collisions, we 
note that their hydrodynamical modeling plays a crucial role in deducing 
the experimental findings [23]. Important sources of uncertainty are given 
by the initial conditions, that are not known precisely. The experimental 
estimation of the initial energy density has to be taken with care, because 
it contains simplifications. In fact, Bjorken estimation of the initial energy 
density at a conservative thermalization time Tq = Ifm/c is used assuming a 
non interacting QGP (i.e. SB equation of state) and a transverse energy dis- 
tribution per unit of rapidity proper time independent during the evolution 
of the system [24, 25, 26, 27]. This second assumption allows to avoid any 
estimation of the lifetime of the deconfined phase but is a strong hj^othesis 
that can lead to changes in the values of the transverse energy density by 
a factor 2 at times between 1 and 8 fm/c [28]. Therefore, already in the 
assumption of a non dissipative QGP the estimation of the initial energy 
density has to be considered with care. The importance of the initial con- 
dition has been shown by the large elliptic flow, V2, in Au+Au and Cu+Cu 
colhsions at RHIC energies [29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40]. 
The centrality, pseudorapidity, and transverse momentum dependences of 
V2 data are described reasonably well by employing the Glauber-type initial 
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conditions and implementing hadronic dissipative effects in ideal hydrody- 
namic models [41]. However, by replacing the initial conditions from the 
Glauber model to the ones expected from a color glass condensate, elliptic 
flow coefficients overshoot the experimental data. This is due to eccentricity 
larger than the ones in the conventional Glauber model. This discrepancy 
strongly suggests effects of viscosity in the QGP [42]. A combined investi- 
gation of both effects, i.e., viscosity/anisotropy effects and the dependence 
on the initial condition is needed to understand the competition between 
these two aspects thoroughly. A measurable quantity that allows such an 
investigation is the entropy density per unit of rapidity dS/dy. In absence 
of viscous effects dS/ dy is a conserved quantity and therefore delivers direct 
information about the initial condition of the system [28], depending on the 
EoS used. On the other hand, if the shear viscosity, ?7s, can not be neglected, 
dS/dy increases, and the details of the dynamical evolution become impor- 
tant. Hydro dynamical calculations use a simplified picture for ?7s and for the 
equation of state. Usually the viscosity to entropy density ratio as tempera- 
ture independent quantity [42, 43, 44, 45, 46] and/or noninteracting SB-EoS 
are used [47, 48, 49, 50]. Only recently, calculations including a schematic 
temperature dependence of rjs/s in the hadronic phase [51] as well as in the 
QGP [52] have been performed. 

Therefore, a consistent description, where a realistic equation of state as 
well as the temperature dependent shear-viscosity derived within the same 
approach is needed to systematically investigate the competition between 
viscosity and initial condition effects on the hadronic yields. Such a consis- 
tent approach is given by the virial expansion, since it allows to derive not 
only a realistic equation of state for the QGP but also to resolve the whole 
temperature dependence of the shear viscosity within a kinetic theory [53]. 
Starting from the results of our model, in this work we will focus on the 
investigation of the role of the viscosity and the initial conditions for the 
evolution and finally for the entropy production in heavy-ion collisions. In 
order to understand the interplay of viscosity, realistic EoS, and initial condi- 
tions we consider a + 1-dimensional time evolution. This simplified scenario 
is appropriate for the early longitudinal expansion stage of the collision. Be- 
cause we focus on the entropy production, which occurs mostly during the 
early stage of the expansion [44], this approximation is completely justified. 
In fact, the authors of Ref. [44] claim, that, for a massless ideal gas, the en- 
tropy production can be calculated to excellent approximation by assuming 
boost-invariant longitudinal expansion without transverse flow during this 
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period. 

The present work is organized as follows: In Section 2 we briefly recall 
the basics of the hydrodynamical equation of motion for a shear- viscous, lon- 
gitudinally boost and transverse translation invariant system. In particular, 
we discuss how the relativistic perfect fluid (Euler) and the Navier-Stokes 
equations of motion can be derived from these equations as a special case. In 
Section 3 we give a brief explanation of the entropy per unit of rapidity in the 
different scenarios. In Section 4 we present our results by discussing in detail 
the role of the viscosity and the initial condition, i.e., formation time, initial 
temperature, and pressure anisotropy for the evolution of the system and 
the entropy production at RHIC at ^/s^^ =130 GeV and the consequences 
for an improvement of their experimental estimation. The conclusions in 
Section 5 finalize this work. 

2. Boost-invariant hydrodynamics 

Relativistic hydrodynamic is based on the local energy-momentum and 
charge conservation 

9,,T^'^(x) = , dXi^)=0, (1) 

expressed in terms of the energy- momentum tensor T^'^{x) and of the local 
charge density n(x). In ideal (Euler) hydrodynamics dissipative effects are 
neglected. Therefore the energy-momentum tensor and the charge density 
are given by 

Tj^ = [e + P)u^'u'' - Pg"'' (2) 
<D = nu>', (3) 

where e and P are the energy density and the pressure, respectively; u'^ is the 
flow four-velocity normalized to m'^m^ = 1 using the standard metric tensor 
giiu _ ciiag(l, — 1, — 1, — 1). The simplest extension to a dissipative regime 
is the introduction of additive corrections linear in flow and temperature 



gradients [47], the so-called Navier-Stokes (NS) approximation 

5T^^ = r/s(V'^n'^ + V^n^ - ^A^'^S^n,) + CA'^'^S^u,, (4) 

- -^(^^)'^'(f). W 

A^'^ = g^^'^-u^'u'' and A^' = A'^'d^, (6) 
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where ( indicates the bulk viscosity, and is the heat conductivity of the 
system. The hydrodynamic equations in NS approximation can be derived 
from a general non-equilibrium theory, the on-shell covariant transport [54, 
55]. The relativistic NS equations are parabolic and therefore acausal. The 
solution formulated by Israel and Stewart (IS) [56, 57] transforms the NS 
equations into relaxation equations for the shear stress, vr^*^, bulk pressure, 
n, and heat flow, q^. The corrections are given by 

71 

ST'"' = tt"" - nA^^ Sn'' = q" (7) 

e + p 

with 

= 0, n^TT^" = u^n"'' = 0. (8) 

The original derivation of the IS equations is not a systematically controlled 
approximation of the transport theory because it is not an expansion in some 
small parameter. In Ref. [57] a quadratic ansatz for the deviation from local 
equilibrium has been employed. Nevertheless, in Ref. [58] a new method for 
deriving the fluid-dynamical equations has been proposed. In this novel ap- 
proach the equation for the dissipative currents are directly obtained from the 
definitions of the currents. We note that, although these equations of motion 
are formally identical to the original IS equations, the coefficients are differ- 
ent. For a detailed overview of the IS theory see Ref. [47]. Here we remark 
that the starting point is given by an entropy current that includes terms 
up to quadratic order in dissipative quantities. These terms are expressed - 
using Landau frame notation - by the same coefficients Oq, ai, /3o, /32, that 
encode additional transport properties. In particular, the set describes 
the relaxation times for dissipative quantities as 

Tn = C/3o , i"q = i^c^TPi , = 2?7s/32 • (9) 

The NS theory is recovered when all these coefficients are set to zero /3o = 

/3i = /32 = tto = tti = 0. 

In the following we focus on a viscous, longitudinally boost-invariant sys- 
tem with transverse translation invariance and vanishing bulk viscosity. As a 
boost-invariant system we consider a system with a longitudinal scaling flow, 
V = (0, 0, z/t), and where all scalar quantities are independent of the spatial 
rapidity, deflned by 

7^ = ^ln[(t + ^)/(t-2)]. (10) 
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If the initial densities are assumed to depend on t and z only through the 
Bjorken (longitudinal) proper time 
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the expansion will evolve such that densities remain independent of 77. The 
Vz will retain the scaling form = z/t [59]. Accordingly, all vector and ten- 
sor quantities can be obtained from their values at 77 = by an appropriate 
Lorentz boost. Because of the symmetries of the system, i.e., longitudinal 
boost invariance, axial symmetry in the transverse plane and rj — > —77 re- 
flection symmetry, the heat flow is zero everywhere. Furthermore, only the 
viscous corrections to the longitudinal, ttl, and transverse pressure, ttt, i.e., 
components of the shear stress tensor evaluated in the 
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vr^z and tt^ 

local rest frame, are non-vanishing. The choice to neglect bulk viscosity is 
sensible since shear viscosity is expected to dominate at RHIC The absence 
of heat flow as well as bulk viscosity leads to 



Tn 



0. 



(12) 



Only the relaxation time, r^r, enters in the equation of motion. With these 
assumptions the equations of motion can be written as 



e + 



It 



n 

h-\ — 
r 


= 


e + p 


_ 


T 




T 


\)] 






VsT J 




3r 


TTt 


2 



(13) 
(14) 

(15) 

(16) 



where the 'dot' denotes d/dr. This special case is a useful approximation to 
the early longitudinal expansion stage of a heavy ion collision for observables 
near midrapidity 77 [47]. Obviously, Eq.(13) describes particle conserva- 
tion and can be solved as 

7-0 n{To) 



n{T) 



(17) 



In Ref. [47] these equation of motions have been studied with the assumption 
of an ideal massless equation of state. Because of this oversimplification the 
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density equation (13) decouples entirely leading to two coupled equations 
for the equilibrium pressure and the viscous correction, vtl- For two limiting 
scenarios, - for the scale invariant case of a constant shear viscosity to entropy 
density ratio, r/s/s, and for dynamics driven by a constant cross section - 
analytic approximate solutions can be found. However, we focus on a realistic 
equation of state and therefore the limitation of a massless ideal system 
must be neglected. We then use the results of the virial expansion for the 
equation of state [22] as well as for the shear viscosity [53], where, within a 
kinetic approach in the ultra relativistic regime, the transport cross section 
as function of the temperature has been consistently evaluated. Therefore, 
we use the relaxation time with the ultra relativistic assumption. As in 
Ref. [57] we obtain 

= -. (18) 

Note that in this way we can systematically investigate the interplay of differ- 
ent effects - realistic equation of state, dissipation, dependence on the initial 
condition- in the same framework. In the following we briefly discuss these 
different scenarios in detail and give the analytic solution of the equation of 
motion, where possible. 

2. 1 . Ideal Hydrodynamics 

The simplest case is that the system does not suffer any dissipative effects. 
Accordingly, we have 

= = 0, (19) 
and the relevant equation of motion is 

s+'-±^ = 0. (20) 
r 

This automatically follows from ttl = vtt = 0. Equation (20) is equivalent to 

^isr) = 0, (21) 

which indicates that st is a conserved quantity. This conservation law is 
broken if dissipative effects emerge. Thus the entropy density evolution is 
given by 

s(r) = with So = s{tq). (22) 
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Consequently, using an EoS a solution of the equation of motion (EoM) can 
be given without explicitly solving Eq.(20). A formal solution for the energy 
density evolution e = e{T) has the form 



— j , with Cq = - and Sq = e{To). (23) 

This results from the separable structure of the EoM, or, more precisely, 
because it is an exact differential equation. Obviously, this formal solution 
holds also for the Stefan-Boltzmann limit, where Cq^ = 1/3 coinciding with 
the sound velocity defined by 

2 dP del P 

This leads to the well known power law £^^(r) oc (tq/t)"^/^. In general for a 
realistic EoS an explicit solution, £:(r), is not possible, because the Cq itself is 
a function of e. However, for the QGP, the ratio pressure to energy density 
can be parameterized using the phenomenological ansatz [60] 

1.\(C^^\, (25, 



3 V l + Be^ 

that provides a good fit to the lattice data in the interval 1.3 < e^^^ / {GeY /iw?Y^^ < 
6 with C = 0.964(5), A = 1.16(6) and B = 0.26(3) fmVGeV [8]. Using this 
parametrization we can explicitly solve the equation of motion (20) as 

F{e) = Fo— with Fq = F{eo). (26) 

T 

The solving polynomial function is given by 

F{e) = e^{e + Lo)^^^ (27) 

with the constants, 

D = j + l-A (28) 
3A D 

naturally, setting A = and C = 1 - and D = 4/3 and u = 1 - the evolution 
in the SB-limit is recovered. 
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2.2. Navier-Stokes 

As mentioned previously, the simplest way to introduce dissipative effects 
is considering the linear Navier-Stokes equations. In this approximation all 
relaxation times are vanishing and consequently the equations for the energy 
density evolution and the viscous correction, ttl, completely decouple 



T^L = (31) 

The NS equation of motion (30) can be reduced to an exact differential 
equation. Formally, using the same notation as before, we can write 

e{T)=so[^y^''°H^s{Vs), (32) 
where the dissipative correction, i^NS? is given by 

HmiVs) = exp ^ dr'^-^l , (33) 

Eq. (33) corresponds to the integration factor needed to transform the NS 
equation motion in to an exact differential equation. This structure can be 
found as well in Ref. [61], where, by assuming for the evolution of the viscosity 

r]f = CnseT + Cnie/T, (34) 

non local corrections governed by the parameter, C^i, have been included. 
With this analytic form for the viscosity the dissipative correction given by 
Htsss can be easily calculated, and one finds 

(l+cg-(4/3)CNs) 



^(^o)(^)"^^^-'""^"'exp ,35) 



T 



as shown in [61] for massless non-interacting particles, i.e., setting Cg = 1/3. 
This assumption for the viscosity means that an additional coupled equation 
for the time evolution of ?7s has been implicitly added in the investigation. 
This is not our strategy, because we have completely determined the shear 
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viscosity in a systematic way, and we can also calculate the time evolution 
without further assumptions. In fact, the EoS and the knowledge of the 
dependence of the viscosity on the thermodynamics quantities,. i.e., on the 
energy density, e [22, 53], allows us to find a coupled system of equations 
for the NS scenario as well as in the general case. Therefore, by finding a 
parametrization of ?7s = r]s{e), we are able to give an analytic solution cal- 
culating H^^lrjs). We note in passing that in this way the equation for the 
energy density and for the viscous correction for the pressure decouple. How- 
ever, we prefer to directly solve the equations of motion without introducing 
any parametrization for the viscosity. 



3. Entropy production 

Now we define the basic observable investigated in this study, and consider 
its different evolution in the different scenarios. In the rest of the paper, the 
subscript '0' refers to the value of quantities at the initial time, tq. 

The most natural quantitative measure of dissipative effects is the entropy 
production. An often used quantity is the entropy per unit rapidity 

^ d^XT s = tAts, (36) 



dy 

where At is the transverse area of the system. This quantity highlights 
several aspects of the QGP because the final entropy production can be 
measured [28] and, additionally, it is sensitive to the dissipative properties of 
the system. In ideal hydrodynamics, the entropy per unit rapidity is evidently 
a constant of the evolution. By the experimental measurement of the entropy 
production at the final time Tf one can conclude the initial condition as 

= ~7~(^o) = tqAtSq. (37) 
dy dy 

In the NS scenario and in the IS regime, dissipative effects determine an 
enhancement of the entropy production. This viscous correction depends of 
course on the initial conditions, specifically on the formation time and initial 
energy density as usual, but in general also on the initial pressure anisotropy. 
Therefore, the most natural choice for the initial conditions are the initial 
entropy density, Sq, and the initial ratio of viscous longitudinal shear and 
equilibrium pressure, defined as 

eo - (38) 
Po 
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A useful equivalent measure is the pressure anisotropy coefficient 



flo ^ f!M = , (39) 

Pto 1 - ^0/2 

which is the ratio of the transverse and longitudinal pressures ptq = p— vrLo/2, 
Plq = Pq + ttlo- Evidently, = indicates central collisions, where pressure 
anisotropy can not develop, whereas peripheral collisions lead to larger |^o|- 
We have shown in Section 2.1 that in ideal hydrodynamics the viscous 
corrections to the pressure always vanish. In this scenario the anisotropy is 
unity and 

= 1 and C = 0- (40) 

are automatically fulffiled. 

As explained in Section 2.2, in the NS approximation the evolution of the 
energy density and of the viscosity are uniquely determinated by Eq.(30). 
Therefore, the initial anisotropy is automatically fixed by the initial condition 
over the entropy density sq and any dependence on (or Ro) vanishes. In 
IS the dependence on the initial pressure anisotropy emerges and thus in this 
scenario we need it as a second initial condition. Additionally, Eq.(36) shows 
that the entropy production is proper time dependent. Therefore, the final 
proper time Tf enters as an important parameter to describe the experimental 
data. 



4. Application to 130 GeV Au + Au collisions 

As mentioned in the previous section the entropy production allows - in 
principle - to quantify the dissipative effects. Furthermore, experimental re- 
sults are available for Au + Au collision at RHIC at a/snn = 130 GeV [28], 
where the final entropy production per unit rapidity is 4501 with an un- 
certainty of about 10%. In the following we attempt to categorize which 
(dissipative) effects and initial conditions are compatible with these experi- 
mental findings. As in Ref. [28] we compare the experimental value to the 
hydrodynamic dS/dy as a function of the initial energy density, Eq- We em- 
phasize that ^0 is strongly dependent on the choice of the formation time, Tq. 
In Ref. [28] the authors use a conservative formation time, Tq = 1 fm/c. 

Our aim is the investigation of the role of all initial conditions and there- 
fore the explicit dependence on the formation time has to be considered. As 
a first step, we focus on the importance to use a realistic EoS, in order to 
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Figure 1: (Color online) The entropy per unit rapidity is displayed as a function of the 
initial energy density eo for the QGP using a realistic EoS [22] (solid lines) in comparison 
to the SB limit (dashed lines) for different formation times, to = 0.1, 0.3, 0.6, 1 fm/c (red, 
green, blue and orange lines) . The horizontal band shows the final-state entropy extracted 
from experiment. Hydrodynamic simulations typically have energy densities between 10 
and 15 GcV/fm'^ and formation times of about 0.6 fm/c. 

clarify the amount and the behavior of the effects that a reahstic EoS gener- 
ates in the experimental observables. In Fig. 1 the entropy per unit rapidity 
is shown as a function of the initial energy density, Eq, for the QGP using the 
realistic EoS of [22] (solid lines ) for various formation times. In following we 
refer to such calculation as non dissipative scenario (ND). For comparison, 
we indicate by the dashed line the SB-limit for the same initialization times 
using the same color coding. The horizontal band shows the final-state en- 
tropy extracted from experiment. Evidently, too fast thermalization times, 
To = 0.1 fm/c and Tq = 0.3 fm/c, lead to an unreasonably high initial en- 
ergy density. Hydrodynamic simulations typically use formation times about 
0.6 fm/c, which corresponds to the blue lines of Fig. 1. The corresponding 
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Figure 2: (Color online) The entropy per unit rapidity is picted as a function of the initial 
energy density, £o, for a dissipative realistic QGP [22, 53] within the NS approximation 
(solid lines) in comparison to the ND results (dashed lines) and the SB limit (pointed line) 
for different formation times, tq = 0.1, 0.3, 0.6, 1 fm/c (red, green, blue and orange lines). 
The horizontal band shows the final-state entropy extracted from experiment. 



results for the entropy in our calculation suggests energy densities between 
10 and 15 GeV/fni'^, which is in agreement with ideal hydrodynamic simu- 
lations. We remark here that for this typical hydrodynamic formation time 
as well as for the conservative one, Tq = 1 fm/c, used in several evaluations 
(extrapolations) of measured experimental data, the difference between SB 
and realistic EoS is not negligible. Therefore, the precise determination of 
the initial conditions is necessary for a good description of the evolution of 
the system. In this context, an implementation of the viscosity effect have a 
larger impact. 

Because of the dissipative effects in the QGP, the entropy density per 
unit rapidity is an increasing function of the proper time. Following [44] we 
assume that the final entropy is mostly produced in the early phases of the 
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collisions, and then we can neglect the contribution of the hadronic phase. 
Therefore the final time, Tf, is automatically fixed, namely as the time when 
the energy density of the system is equal to the critical energy density. 



Because in the NS approximation the initial viscosity of the system is com- 
pletely determined by the initial energy density, we can directly compare the 
results for the entropy per unit rapidity for the non dissipative system shown 
before with the ones from the NS equation of motions. Using the same color 
coding as before, we show in Fig. 2 the entropy per unit rapidity as a function 
of the initial energy density Eq for the dissipative QGP [22, 53] within the NS 
approximation (solid lines) for various formation times in comparison with 
the non dissipative and SB calculations. For all formation times one observes 
an enhancement of the entropy production by the NS viscous corrections in 
comparison to the ND results. For tq = 1 fm/c and tq = 0.6 fm/c we find 
moderate corrections, that seem to recover the results of the SB ideal fluid 
shown in Fig. 1. In particular, for tq = 0.6 fm/c, the effect of a realistic 
equation of state and the linear (NS) viscous terms compensate also quanti- 
tatively each other. For the faster thermalization scenarios, r = 0.3 fm/c and 
r = 0.1 fm/c, the enhancement is very significant and allows to describe the 
experimental finding for the entropy per unit of rapidity using smaller - but 
nevertheless high - Eq of about 20 GeV/fm^ and 25 GeV/fm^ respectively. 
The formation time, Tq, indicates an atypical behavior, not only because of 
the pronounced enhancement of the entropy production, but primarily be- 
cause of the change from a concave to a convex function in dS/dy{Eo)- This 
can be a hint that results in NS approximation with small formation time 
have to be questioned. Deciding this question, one has to consider the full 
IS equations of motion for the entropy production. 

As explained in Section 3 the initial conditions for the IS equation of 
motion are not completely fixed by the initial energy density of the system 
because of the additional explicit anisotropy of the system. Therefore, one 
has a family of solutions labeled by the anisotropy coefficient, ^q, as additional 
parameter. For each formation time, we calculate, for different values of the 
initial anisotropy, the entropy production at T = Tc as function of the initial 
energy density. We use i^o = 0, = —0.5, and = —1, which correspond 
to the values for the anisotropy coefficient of Rq = 1, Rq = 2/3, and Rq = 0, 
respectively. We prefer to concentrate on two limiting cases, tq = 0.1 fm/c 
and Tq = 1 fm/c. 




(41) 
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Figure 3: (Color online) The entropy per unit rapidity is depicted as a function of the 
initial energy density, Eq, for a dissipative realistic QGP [22, 53] at tq 0.1 fm/c within the 
IS approach for different values of the initial pressure anisotropy, ^o = 0,— 1,— 2,— 3 (black, 
turquoise, violet and brown solid line respectively) in comparison to the NS approximation 
(red solid line), the ND scenario (dashed red line) and the SB limit (dotted red line). The 
horizontal band shows the final-state entropy extracted from experiment. 



In Fig. 3 we show the entropy per unit rapidity as a function of the initial 
energy density Eq with the formation time Tq = 0.1 fm/c within the IS ap- 
proach. The different values of the initial pressure anisotropy = 0, —0.5, — 1 
are displayed as black, violet and turquoise lines, respectively. For compar- 
ison we added the corresponding results within the NS approximation (red 
solid line), the ND scenario (dashed red line) and the SB limit (dotted red 
line). The experimental result is located within the horizontal band in Fig. 3. 
Already at vanishing anisotropy the deviation from the previous results is ev- 
ident in the whole energy range, the difference to the non dissipative scenario 
are small for low initial energy density, e > 10 GeV/fm^. By increasing the 
anisotropy of the system, the entropy production also increases. For all val- 
ues of .^0 we observe a linear relation between dS/dy and Eq. Nevertheless, we 
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Figure 4: (Color online) The entropy per unit rapidity is displayed as a function of the 
initial energy density, eo, for a dissipative realistic QGP [22, 53] at tq = 1 fm/c within the 
IS approach for different values of the initial pressure anisotropy, = 0, —1, —2, —3 (black, 
turquoise, violet and brown solid line respectively) in comparison to the NS approximation 
(orange solid line) and the ND scenario (dashed orange line) and the SB limit (dotted 
orange line). The horizontal band shows the final-state entropy extracted from experiment. 

reject such a fast formation time for two reasons: not only does the NS formu- 
lation not allow To = 0.1 fm/c, but also at extremely high initial anisotropy 
^0 = —1, that corresponds to the limiting case with the tranverse pressure 
only, high values of the initial energy density are needed to reproduce the 
experimental entropy production. That makes this scenario unlikely. The 
contrary situation is given by a conservative formation time, tq = 1 fm/c. 

In this case Fig. 4 shows the results for the entropy per unit rapidity as a 
function of the initial energy density Eq within the IS approach for the same 
values of as before. Again, we compare to the corresponding results in 
the NS approximation (orange solid line), the ND scenario (dashed orange 
line) and the SB limit (dotted orange line). We note the same quahtative be- 
havior as for the small formation time (Fig. 3): by increasing the anisotropy 
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Table 1: Results for the deviation from the SB estimation of the initial energy density in 
the ND and NS scenario 



To 

[f m / c] 


^SB 

[GeV/fm3] 


^ND /^SB 
^0 /^O 


c-NS /c-SB 
^0 /^O 


0.1 


143.10 


1.05 


0.20 


0.3 


33.07 


1.10 


0.73 


0.6 


13.12 


1.15 


0.99 


1.0 


6.64 


1.21 


1.14 



of the system, the entropy production also increases. In particular, we find 
an initial energy density between Eq ^ 7 — 9 GeV/fm"^, that is larger than 
experimental estimation in Ref. [62] of about eo ~ 4 — 5 GeV/fm^. This 
is not surprising, because the experimental measure is based on very sim- 
plified and strong assumptions, that do not give a realistic description of 
the system. This suggests that the estimation procedures have to be im- 
proved because of the effects of the final proper time Tf, the viscosity and 
the anisotropy can not be neglected. To be more quantitative we can ex- 
tract from the experimental estimation of the final entropy the initial energy 
density for different formation times in the different scenarios. In this way 
we can compare the (typical experimental) values, ^q^, obtained using the 
ideal Bjorken expansion, i.e., the SB limit, with the values Sq^, Eq^ and e^q 
in the ND scenario, in the NS and IS theory respectively. Evidently, for the 
IS calculation, non only the parametric dependence on the formation time 
tq, but also on the initial pressure anisotropy has been considered. The 
results are summarized in Tab. 1 and Tab 2, where the deviation from the 
standard SB estimation of the ND, NS and IS calculation is listed. For small 
formation times (tq = 0.1 fm/c and tq = 0.3) the effects of the relativistic 
equations are almost negligible, i.e. lower then 10%. However, dissipative 
effects are to be included, because the viscosity leads to an deviation be- 
tween 30% and 80% in the NS scenario and between 40% and 60% in the IS 
scenario. Only for central collisions (^o = 0) the deviation in the formation 
time Tq = 0.3 fm/c is maybe reasonable (10%), although, in this case, the 
assumption leads to an unlikely large initial energy density. Therefore, such 
small formations time have to be rejected. For tq = 0.6 fm/c, often used 
in hydrodynamical calculations, we note that the viscous effects completely 
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Tabic 2: Results for deviation from the SB estimation of the initial energy density in the 
IS scenario 



To 

[fm/c] 


^SB 

[GeV/fm3] 


e-IS /c-SB 
^0 /^O ' 

^0 = 


^IS /^SB 
^0 /^O ' 

^0 = -0.5 


^IS /^SB 
^0 /^O ' 

eo = -l 


0.1 


143.10 


0.62 


0.50 


0.41 


0.33 


3.07 


0.9 


0.80 


0.72 


0.6 


13.12 


1.02 


0.99 


0.93 


1.0 


6.64 


1.16 


1.12 


1.08 



compensate the correction of the reahstic EoS in the NS scenario and in the 
IS theory for non peripheral colhsions. Thus, in this case, using the Bjorken 
expansion can be justified. For the conservative formation time, tq = 1 fm/c, 
the effects of the reahstic EoS are more important, about 20%, and can not 
be compensated by the inclusion of NS viscous corrections. Therefore, it is 
not surprising, that also in central collisions in the IS scenario the deviation 
is also sizable. For more peripheral collisions the SB limit seems to lead to a 
more or less satisfactory approximation (about 8%) for the IS calculations. 
Nevertheless, from this discussion it is evident, that the use of the Bjorken 
expansion for the evaluation of the initial energy density, without consider- 
ations of interplay between formation time, viscosity, and the centrality of 
the collisions can be very questionable. Additionally, the quite satisfactory 
agreement between IS hydrodynamics calculations and the transport ones 
within the cascade BAMPS for small ?7s [63, 64] also remarks the validity of 
the IS scenario. 

5. Conclusions 

In this work we have discussed the importance of the initial condition and 
the role of the viscosity for the evolution of the fireball of heavy-ion collisions. 
We have used a realistic EoS derived within a virial expansion [22], that is 
in line with recent three- flavor lattice QCD data [8]. The shear viscosity has 
been consistently calculated within this formalism using a kinetic approach 
in the ultra-relativistic regime with an explicit and systematic evaluation of 
the transport cross section as a function of temperature [53]. We explicitly 
considered different scenarios: ideal hydrodynamic, dissipative effects in the 
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Navier-Stokes as well as in the Israel-Stewart formalism, from conservative to 
very fast equilibration dynamics. We choose the parameter of these studies 
in order to describe the experimental findings of the entropy production for 
Au + Au collision at RHIC at a/snn = 130 GeV. The assumption of a fast 
equilibration would require unreasonably high values of the initial energy 
density. With the conservative (tq = 1 fm/c) and the typical hydrodynam- 
ical (to = 0.6 fm/c) formation times the initial energy density needed to 
reproduce the final entropy is more in line with the experimental estima- 
tions. In these scenarios, the interplay between effects of the viscosity and 
of the realistic EoS can not be neglected in the reconstruction of the initial 
state from the experimental data. Additionally, centrality dependence have 
to be considered, because different impact parameters lead to different initial 
anisotropy values, which modify sizeably the estimation of the initial energy 
density. In conclusion, our investigation shows that from the experimental 
final entropy it is very hard to derive unambiguous information about the 
initial conditions and/or the evolution of the system. The choice of the for- 
mation time, the viscosity, and the initial anisotropy are interplaying aspects 
which have to be included in the estimation of the initial condition. There- 
fore, we suggest that the easy extraction rule used should be revised for a 
a better estimation of the uncertainty in the measurement. Of course, this 
improvement is model dependent, but we note that also the usual evaluation 
using the Bjorken expansion is a model calculation. 

Furthermore, the solution of the IS equation with a realistic equation 
of state and with the inclusion of viscosity leads to further interesting ap- 
plications in the description of the heavy-ion collisions, as (semi) analytic 
solution of the hydrodynamical equation in accord of Ref . [65] . Additionally, 
a systematic calculation of the bulk viscosity within the virial expansion and 
the kinetic theory can be implemented to consider all remaining dissipative 
effects in the QGP in a systematic way. 
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